HBT correlation in 158 A GeV Pb+Pb collisions 
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The large acceptance TPCs of the NA49 spectrometer allow for a systematic mul- 
tidimensional study of two-particle correlations in different part of phase space. 
Results from Bertsch-Pratt and Yano-Koonin-Podgoretskii parametrizations are 
presented differentially in transverse pair momentum and pair rapidity. These stud- 
ies give an insight into the dynamical space-time evolution of rclativistic Pb+Pb 
collisions, which is dominated by longitudinal expansion. 



1 Introduction 

Recent high statistics experiments have demonstrated that, in heavy ion colli- 
sions at relativistic beam energies, a single characteristic radius from intensity 
interferometry does not exist. First of all one deals with systems, in which 
the emission of particles is distributed over various ranges in three dimen- 
sional space as well as in time. Secondly, due to the dynamical behaviour of 
the source (eg. collective expansion) correlations arise between the momenta of 
particles and the point of emission in space-time, which results in a dependence 
of HBT radii on kinematicai quantities. In the case of a "boost invariant" pic- 
ture, first suggested by Bjorkerfl, the hot and dense initial stage of the collision 
expands longitudinally until freeze-out such, that a distinct longitudinal ve- 
locity (v z ) profile is established: v z — z/tf. Here, z is the distance of particle 
emission at freeze-out from the center of the collision and tf the freeze-out 
time. With respect to HBT correlations - appearing at low momentum differ- 
ences of particle pairs - such a velocity profile disconnects different parts of 
the source along the beam axis. Only additional thermal motion with average 
velocities of < Vtherm >~ y/T/m± can - over certain distances - compensate 
for the longitudinal velocity gradient. This means that, in case of relativistic 
heavy ion collisions, HBT radii do not measure the geometry of the source, but 
a "length of homogeneity" □ dz = R z — tf \jT/mi_, which may vary inside the 
source and hence depend on the phase space in which the pairs are observed. 

Experimentally, the question of space-time size of the source is addressed 
by extracting HBT radii for all space-time components of the momentum dif- 
ference vector q 11 = p 1 ^ — pt^in = 0..3 of pairs of identical particles 1 and 2. 
Unfortunately, the "on-mass-shell"- constraint for the observed (real) parti- 
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cles reduces the degrees of freedom to three. For pairs this constraint can be 
written as 

q k = qk (1) 

with the average pair momentum fc M = (p^ +p%)/2. Later it will be shown how 
this condition is utilized to eliminate one of the four space-time components 
of q^ L to deduce parametrizations of the three dimensional correlation function 
^jM?) according to Bertsch-Pratt(BPu) or to Yano-Koonin-Podgorctskii (YKP 

The second point - the dynamical space momentum correlation - is inves- 
tigated experimentally by carrying out full three dimensional HBT analyses 
for intervals in the average transverse pair momentum and average pair 
rapidity (y — (yi + y2)/2), separately. All together this can be expressed by 
(q components introduced in section |3|) : 

Ca(«") — = C 2 

Such a detailed analysis requires of course a large data sample over a large 
fraction of phase space. Facilitated by the large number of particles produced 
in central Pb+Pb collisions at 158 A-GeV and the large acceptance of the 
NA49 spectrometer at the CERN SPS accelerator such a study has become 
feasible and is presented here. 

2 NA49 spectrometer and data analysis 

The NA49 spectrometeii (figure [l]) is located in the North area of the CERN 
SPS. The results presented here are derived from the run period in 1995 with 
a 158 A-GeV Pb-beam. In every spill (5 s) about 10 5 Pb 82+ ions impinge on a 
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224 mg/cm 2 Pb-target. 1% of these projectiles interact with the nuclei of the 
target foil. The most central 5% (impact parameter b = — 3.5 fm) of these 
collision are selected by the NA49 trigger utilizing a Veto- Calorimeter, which 
measures the energy carried by the beam spectator fragments. In such central 
collision about 1200 charged particles are produced, roughly 800 of them are 
detected in at least one of the four TPCs (Time Projection Chambers) of 
NA49. At 2.0 m and 5.8 m downstream of the target two smaller chambers 
(VTPC1, VTPC2; volume=2.0 x 0.7 x 2.5 m 3 ) are placed inside magnetic 
dipole fields. In the "Standard Configuration" the fields amount for 1.5 Telsa 
(VTPC1) and 1.1 Tesla (VTPC2), respectively, whereas in case of the "Low 
field configuration" , fields of 0.3 Tesla and 1.5 Tesla have been used, shifting the 
acceptance for pions towards central rapidities. Further downstream - outside 
the magnetic field - the acceptance is extended to forward rapidities and higher 
momenta by two large volume TPCs' (MTPC; volume=3.8 x 1.2 x 3.8 m 3 ). 

All four TPCs are equipped with a charge sensitive read-out, highly seg- 
mented along the beam (z-axis) and perpendicular to it in horizontal direction 
(x-axis). For each of these 180.000 pads the charge is determined in 512 consec- 
utive time slices (100 ns), to determine the vertical (y-axis) position of particle 
tracks via the drift velocity (v^ ral.4 cm//is inside VTPC and w2.2 cm//xs in- 
side MTPC). This scheme allows for 3-dimensional track reconstruction. From 
the curvature of reconstructed tracks inside the VTPC the momentum of par- 
ticles is determined to a precision of Sp/p 2 ps 0.3% in addition to their charge. 
Moreover, all chambers measure the energy (dE/dx) deposited by the particle 
in the detector gas, which - in a future stage of analysis - will be used to 
identify the particles. For the HBT analysis presented here, pairs of hadrons 
of the same charge (h~h~ or h + h + ) without further identification have been 
considered. These are dominated by pairs of identical pions (7r _ 7r~ or tt + tt + ). 
It has been shown, that a contamination from other particles has negligible in- 
fluence on the determination of HBT radii, with the exception of the chaoticy 
parameter A, which will not be discussed in this article. 

Independent analyses have been carried out for the VTPC2 and the MTPCs. 
The VTPC2 analysis is based on 40.000 central Pb+Pb collisions, half of them 
taken in the "Low field configuration" and the other half in "Standard field 
configuration". The analysis of the MTPC data] includes the same 40.000 
events and adds another 50.000 events in "Standard field configuration". In 
both analyses pairs of tracks with distances less than 2 cm inside the TPC 
have been excluded to eliminate the influence of two particle reconstruction 
inefficiency for very close tracks. The same requirement has been imposed on 
the distribution of uncorrelated pairs, which is generated by combining tracks 
from different events. This distribution is used as reference (denominator) in 
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the determination of the correlation function. Both data sets have been cor- 
rected for Coulomb final state interaction based on measured correlations of 
opposite-sign charged particles. In case of the VTPC2-d&t& this correction is 
carried out by the Gamov factor G(qi nv ); qi nv = \/—(jpx pr p2) 2 modified by a 
dumping term to account for the finite size of the sourceu: 

C 2 — rr = C^- eas , x [(G(q mv ) - 1) C -«W*« + l) (3) 

The parameter q e ff is determined from a fit of the 1-dimensional correlation 
function C\ (qinv)- 

The MTPC-data are corrected by the correlation function C\ of the 
opposite charged particles, evaluated in the same three dimensionial relative 
momentum space as the like-sign pairs; eg. for BP-paramcrtization (see below): 

^2.corr = ^2,meas. X ^ + (iside, Qout, Qlong) (4) 

No further corrections have been applied to the data. The systematic 
error in the determination of the HBT radii is estimated to about 7%. The 
correlations function have been obtained in different reference frames. The 
most intuitive choice of a reference frame for all pair momenta might be the 
Center- of Mass System (CMS) of the colliding ions, but in case of variations 
of the longitudinal velocity across the source due to a longitudinal expansion, 
the Longitudinal Co-Moving System (LCMS) might be better suited. It is 
defined on a pair-by-pair basis such that the longitudinal pair momentum fcii 
vanishes. The Fixed Longitudinal Co-Moving System (FLCMS), in which 
the observer frame for every interval in pair rapidity is fixed at the center of 
that interval, might be seen as a compromise between both. In case of narrow 
widths of rapidity intervals LCMS and FLCMS are equivalent and are therefore 
treated as one in the following discussion, even though the MTPC dataQ have 
been evaluated in the CMS and FLCMS frames, whereas the VTPC2 datatJ 
use the CMS and LCMS frames for BP projections and CMS and FLCMS in 
case of YKP. 

3 Q- parametrizations 

As shown in references 00, the correlation function for a source distribution, 
expanded to second order at the (in general fc-dependent) space time points of 
maximum emission (X; saddle point) can be written as: 

C 2 = 1 + e -9 M 9"<*M*»> (5) 

with = — X^; <> denotes an averaging over the source distribution. 
Even though this "model independent" expression (0) can - because of the 
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Figure 2: 1- and 2-dimensional projections of the h~h~ correlation function measured 
in the N A49 MTPC using BP parameters for pair rapidities 4 < y < 5 and transverse 
momenta k± < 100 MeV/c (FLCMS). The projections on components qj have been 
carried out by restricting the other components qn i / j to qt < 30 MeV/c. 



constraint in cq. (gj) - not be used to describe measured correlation functions, 
it is a generalization of both the Bertsch Pratt B as well as the Yano Koonin 
Podgoretskii parametrization. It therefore gains its importance as a tool when 
interpreting the radii extracted in the framework of both formulations; it fur- 
thermore provides a proof of consistency, when comparing results from both 
parametriztions . 

3. 1 BP parametrization 

Eliminating the temporal component of (|^) by condition (|l|) in the form q = 
flq ((3 —Js/ko), results in a correlation function parametrized according to 
reference H: 

C 2 = 1 + Ae-«»^' ! '^- 9 ™' s »'-«'»» R ?»»- 2 9»'9'.»» s L-,„» s ( 6 ) 

Here qi ong is the component of q 11 in beam direction, whereas q S id e and q ut 
are those perpendicular to it, with q out \\k and q S ide -L k. 

Due to symmetries of the sources considered here, only the "out-long" 
cross term remains, all others ("out-side" and "side-long") vanish. This is 
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consistent with the experimental observation presented in figure || for h~h~ - 
pairs at rapidities 4 < y < 5 and transverse pair momentum kj_ < 100 MeV/c. 
The contours lines in the 2-dimensional projection deviate from a circular 
shape (no cross term correlation) only in case of the "out-long" projection. 
The 1-dimensional projections demonstrate the good agreement between the 
function (|J) and the data for small q, where the correlation signal is clearly 
visible, as well as for larger q, where the data points are consitent with Ci = 1 
within the error-bars. At this point it should be emphasized, that all HBT 
radii presented here are derived by a simultaneous fit of all three components 
of q in Ci. The projections are generated for better visibility only. 

When interpreting the BP radii it is advantageous to explicitly write down 
the relation between equation ^ and ansatz ^[ 

<x\> (7) 

< (x x - (3 J) 2 > (8) 

< (x z -^i) 2 > (9) 

< (x x -[3j)(x z -/?||t) > (10) 

With the exception of R S ide, BP-radii mix spacial (x xyz )a,nd temporal com- 
ponents (t) of the source and an interpretation becomes therefore reference 
frame and model dependent. In case of a longitudinally expanding source 
the "length of homogeneity" observed in the center of mass frame, appears 
Lorentz contracted in different intervals of pair rapidity. Such a behaviour is 
supported by the rapidity dependence of the measured radius Ri ong , as shown 
for different intervals in k± in figure |^. The data points can be described by 
Riong — t f /cosh(y) y/T/m± , which is the Lorentz frame dependent expression 
of section 1. Assuming a temperature of T = 150 MeV, a freeze-out time of 
tf ~ 9 — 7 fm/c can be derived for the different fc^-intervals. By comparing 
the measured "side"- and "out" -radii to equation (j|)-(0) the duration time of 
freeze-out yields m 2 — 4 fm/c. Both radii ("side" and "out") appear to be con- 
stant over y < 2.5 within errors. In addition to the physical interpretation it 
is important to note the good overall agreement between the different analyses 
and between pairs of opposite total charge, i.e. h~~h~ compared to h + h + . 

3.2 YKP parametrization 

Instead of eliminating the temporal component of (^J) in the BP formalism, 
condition (|l|) might also be used via the relation q x = qo/(3± — qn/3n/P±. With 
the picture of a boost invariant source in mind with different parts of the 
source moving at different longitudinal velocities, one explicitly introduces a 
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Figure 3: The dependence of BP radii on the pair rapidity for h~h~-pairs (*) and 
ft+/i + -pairs (+) in the VTPC2 and fe~ft~-pairs (v) in the MTPC in four intervals of 
transverse momentum kt (CMS). The rapidity scale is shifted to the center of mass 
system of the ions. The horizontal error bars correspond to the width of the intervals 
chosen in the analysis and the vertical to the statistical errors only. 



8 




0.1 0.2 0.1 0.2 -0.2 0.2 

Q (GeV) Q ± (GeV) CL (GeV) 




Q ± (GeV) Q (GeV) Q (GeV) 

Figure 4: Projections of the same data sample as in figure ^, but this time in the YKP 
parametrization (FLCMS). The projection onto component(s) qj have been carried 
out by restricting the other component(s) qi\ i 7^ j to qi < 70 MeV/c. 



longitudinal velocity parameter (Pykp) into the correlation function of YKP 
typett 

C 2 = 1 + \e~ q2±R2± ~rt K p(qii-f 3 y-Kpqo) 2 R^-'iYKp(<io-f3YKpq\\) 2 Rl Ql) 

with jyk p = 1/ \/l — Pykp ' 9-L = \J + Qy f° r the transverse momentum 
difference, go for the energy difference and gii for the longitudinal component. 
In this case the interpretation of extracted radii becomes more evident, since 
space and time components are decoupled (the validity of the approximation 
is discussed in referenced): 

R\(k) = <x 2 y > (12) 
R 2 (k) m < P > (13) 
R\ {%) « < > (14) 

Figure ^ shows one and 2-dimensional projections of the measured corre- 
lation function in YKP coordinates evaluated in the LCMS frame for the same 
kj_ and y-interval as in figure |^. Again, the data are described well by the cho- 
sen gaussian ansatz. A problem of this type of parametrization manifests itself 
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Figure 5: The dependence of YKP-HBT radii on the transverse pair momentum k± in 
different intervals of pair rapidity (LCMS) (otherwise same conventions as in figure H). 



in the projection of go- For a given interval in q± and q\\ only a limited region 
of qo is kinematically available. This is the reason for the large uncertainties of 
Rq in figure |5|, which sumarizes the /^-dependence for the YKP-radii in differ- 
ent intervals of rapidity. The clear decrease of vs. k t again points to strong 
space-momentum correlations in the source. Moreover, even in the transverse 
direction -fiLi decreases for larger fc t -values, which can be interpreted by trans- 
verse flow E£L The estimate of the duration of emission given in section 3d is 
confirmed by the extracted values of Ro- 

In a parametrization of YKP type one can gain further insight into the 
dynamics of the source by utilizing the YKP velocity Pykp- The YKP rapidity 
Vyk p = \ In i^y^p + V derived from the measured Pyk p is compared to the 
pair rapidity in figure ^. The data show the characteristics of a source, which 
expands in longitudinal direction. Even though deviations from an ideal boost 
invariant picture (line in figure 0), become apparant at forward rapidities, the 
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Figure 6: Dependence of the Yano Koonin rapidity vs. the pair rapidity for four 
intervals in k± (LCMS). 

consitstency which such a model is good. 
4 Conclusion and Outlook 

Results from different analyses of multidimensional Bcrtsch-Pratt and Yano- 
Koonin-Podgoreskii parametrizations of two particle correlation functions in 
158 A-GeV Pb+Pb collisions have been presented differentially in pair rapidity 
y and transverse pair momentum kt ■ The results confirm the importantance of 
extracting HBT-radii seperately in different parts of phase space to disentangle 
the dynamical correlations of the source. The pion source appears to expand in 
a close to boost invariant way, as seen in the rapidity dependence of the Yano- 
Koonin velocity as well as of the Ri ong radius in the CMS system. Moreover, 
a finite duration time of emission of 2 — 4 fm/c and a decreasing tranverse 
radius at large k t are observed. For a better understanding of the behaviour of 
the radii a comparison to a similar analysis of proton-proton and proton-lead 
collisions is currently in progress. Moreover, the centrality dependence and 
beam energy dependence might constrain interpretations even further and will 
be investigated in upcoming analysis and further data taking. 
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